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ABSTRACT 


The main objective of this thesis is to understand tiie role of pkLstic zone at 
crack tip on mixed mode crack propagation, effect of lexuirc in the drawn steel wire 
and the effect of warm pre-stressing on crack toughness Firstly, crack-path instability 
IS studied in a drawn steel wire . Part-through semi-ellipiic crack surface is observed 
while the crack propagates along 45° helical path. The maximum shear stre.ss elevation 
changes its direction as the crack propagates. Secondly, unstable crack propagation 
IS studied through shear fracture along the axis of the ware. Thirdly, the effect of 
anisotropy on crack tip plastic flow is studied on textured 7075-T6 aluminium alloy 
plate shear strain localization and non-homogeneous plastic flow produce a critical 
shear stress concentration on 45° plane. Finally, Warm pre-stressing effect has been 
studied m mild steel. Pre-cracked specimens are loaded in three point bend test up to 
the point of 60% to 95% of its fracture load at room temperature, then it is fractured 
at different sub-zero temperatures. Crack-toughness is seen to increase due to warm 
pre-stressing. 
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Chapter 1 


UNSTABLE CRACK 
PROPAGATION IN DRAWN 
STEEL WIRE, A STUDY ON 
CRACK TIP PLASTIC FLOW 

1.1 SUMMARY 

The mam purpose of this investigation is to understand the role of plastic anisotropy 
on crack propagation, which normally leads to an unstable fracture in textured materials. 
Firstly, crack-path instability is considered, in reference to, the transition for different 
changes in crack propagation modes 0.6 percent carbon steel drawn wire is subjected 
to combined mode of loading rotation in the presence of an axial tension. Mixed mode 
crack-propagation is observed by SEM with respect to the grain orientation direction of 
the wire. Part-through semi-elliptic cracks are seen to form in a wire while the crack 
propagates along 45° helical path. The maximum shear stress elevation within the plastic 
zone continuously changes its direction with respect to the angle of crack propagation. 
Secondly, an abrupt crack extension is seen to occur due to the plastic anisotropy lead- 
ing to an unstable crack propagation through shear fracture along the axis of the wire. 
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Calculations based on Irwin-sih-srnith overestimated the direction of maximum shear 
stress vector by 3 3° with respect to the axis of the wire Finally, the effect of anisotropy 
on ciack tip plastic flow is studied on textured 7075-TG aluminium alloy plate Experi- 
mental observation shows shear stiain localization These obser^atlons favourably agree 
with the previous work 
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1.2 INTRODUCTION 


Steel wires are normally made by drawing process, so that, grains are seen to be 
elongated along the axis of the wire Such steel wires are extensively used in engineering 
components, for example as valve springs, wire ropes in elevators, bicycle spokes, or even 
as musical instruments. In service wires are subjected to severe vibrations, obviously, 
producing fatigue failure. Fatigue cracks are usually seen to start from a defect on the 
wire surface 1 — 3j. Crack propagation mainly depends on the type of loading on the 
wire Final fracture occurs due to an abrupt propagation of a micro-crack along the axial 
direction of the wire [4 — 5]. This is mainly due to sudden loss of resistance to shear 
deformation along the grain orientation direction in a textured material [6]. Thus, there 
are three distinct stages: firstly, the initiation, which is mainly due to a surface defect 
secondly, the crack propagation, which depends on the type of loading, for example, a 
helical crack is seen to form due to a combined mode of loading under tension and torsion, 
and finally, micro-fracture can suddenly initiate within the plastic zone in the presense 
of anisotropy, (and not on the yielded-unyielded boundary, as predicted by elasticity 
theory [7—11;). It is observed that, the gram flow direction plays an important role in 
micro-fracture imtialisation in metals. 

This project is divided into three parts. The first part deals with the theoritical 
derivation on mixed mode loading on an inclined elastic crack in a circular wire. An 
explicit expression for stress intensity factor is given for a semi-elliptic part- through 
crack on 45° inclined plane [12 — 13]. Slip-line field solutions are then employed to 
calculate the stresses within the yielded zone. The second part is on the study of 
fatigue-fracture transition in 0 6 carbon cold drawn steel wire ( 2mm ) . The last part is 
mainly on crack tip plastic flow in heavily cold rolled 7075-T6 aluminium alloy (4 mm 
thin plate). 

Under the action of combined mode of loading in tension-torsion, fatigue cracks were 
grown in wires from notch Fig. 1(a). Scanning electron microscope was found to be 
suitable for measuring the propagation of helical crack in steel wire (see Fig. 2(a)). Fig. 
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1(b) is a concept of plastic anisotropy m a flat plate under-tension. Fig. 2(b) shows a 
pin loaded aluminium specimen in tension The tip was opened by a static load after 
fatigue cracking Thin film carbon replica micro-grids were glued onto the specimen 
surface to measure the shear strain distribution within ;isymmetric plastic enclave of the 
crack tip. Thus, the objective of the present work is to find the reason for crack path 
instability in anisotropic drawn steel wires, in reference to the mode transition in crack 
propagation. It is observed that, at first the crack propagates in pure mode I then in a 
combined mode of mode I and mode II along the 45° plane, and finally, an abrupt mode 
II pure shear fracture takes place along the axis of the wire. Experiments on (7075-T6 
heavily cold rolled) aluminium alloy gives an idea on the nature of plastic flow at crack 
tip in anisotropic material, but this may not give an exact stress field for a drawn steel 
wire. Hence, for an inclined crack in a wire, the reason for this sudden change in crack 
propagation behavior will be investigated. 

\ erbaak was the first to emphasize that there is a need to make such a careful study 
to prove the conditions favourable for an abrupt change in crack-propagation path. It 
is well known that an inclined crack can suddenly bifurcate under a combined mode of 
loading, even in an ideal isotropic elastic material. This is due to the fact that, the 
direction of the maximum shear stress at the tip of the crack is continuously changing as 
the crack propagates. Finally, the energy minimum condition has to be satisfied for an 
instability Thus, there are three distinct aspects which need to be considered carefully: 
firstly, due to anisotropic texture, the location of the maximum shear stress elevation 
within the plastic zone changes as the crack propagates, secondly, due to the combined 
mode of loading, the direction of the maximum shear stress also changing -nfith respect 
to the angle of crack orientation, and hence finally, at some stage of propagation, the 
direction as well as magnitude must coincide with the lowest resistance offered by the 
wire material to have an instability. 
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1.3 INCLINED PART-THROUGH SEMI-ELLIPTIC 
CRACK IN A WIRE UNDER MIXED 
MODE LOADING 


One rarely finds that an initial starting crack is through the cross-section of a shaft. 
Figs 3(a) and 3(b) show a typical shape of a starting crack in a circular wire For 
the basic mathematical formulation purpose, a simplified approach has been followed as 
suggested by \erpoest, Aernoudt et al [5]. A semi-elliptic part-through crack is assumed 
to be situated on 45° inclined plane in a circular shaft. The co-ordinate systems are 
explained in Figs. 3(a) and 3(c). Irwin [7], Newnan et al [12], and Smith et al [13] have 
analyzed this kind of problem in a rectangular beam. The stress intensity factors in 
Mode I and Mode II will be determined in this paper. The SIFs calculated are valid at 
the tip of the ellipse. 

In cylindrical co-ordinate, for asymmetric situation, the equations of equilibrium, in 
terms of the displacements are [17 — 18] 


V^u + 


Vh + 


1 dA 

\ — dr 
1 dA 

l-2u Vm 


1 .2 dv u 

r T do r 

1 .V 2 du. 

T r r do 

r-,2 1 dA 

V^w + — 

1 — 2u dz 


0 

0 

0 


where, u=radiaJ displacement, v=tangential component of the displacement and w = 
axial displacement along z direction. denotes the Laplacian operator and u is the 
Poisson’s ratio. A being the dilatation, which is given by 

A du u I dw 

^ ~ dr ' r ^ r de dz 

The equations of equilibrium are satisfied, if one takes 


u = 

V = 

w = 


d'^cj) 2 d-ip 

drdz r do 

d'^(f) 2d'ip 


rdOdz dr 
2(l-p)VV- 


d‘^4> 
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where, = 0, = 0 

Thus, o IS bi-harmonic and ip is harmonic 

The stress components can be written cis 
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where, G = Shear modulus of the wire material. 

One can write the biharmonic function p and the harmonic function ip in the following 
forms. 


p{r,9,z) 

^m=0,l,2, 

[(pmir, z) cos 

m9 — 

cpm{r,z) sin m9 

<P{r, 9, z) 

II 

M 

0 

1 — • 

(pm{r, z) cos 

m9 


ipir,9,z) 

Y'co 

^m=0,l,2, 

[i’mir, z) sin 

m9 — 

i’miT-, z) cos m^j 

ip{r,9,z) 

Y^oo 

ipm{r, z) sin 

mO 



This IS due to the fact that, considering only single valued of m, without loss of generality, 
one can always write (pm = ipm = 0 

It is easily seen that (pm and ipm are solutions of the following partial differential 
equations 



4^m 
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^ m 

9r2 
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dr 


+ 
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dr 
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Using Bessel function from snedden [18] 


r J„Kr) dr 


■ii 

^ dz2 


/o“ Jmi^r) dr = 0 
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for 0 < r < 1 

and 

/» r Vi,h„ J,„(?r) <ir = (i - /“ r cv J.„(s'r) dr = 0 

for 0 < r < 1 

If one transforms \v r.t. I-miQr), then put Q = and 
t'^Im{C,r) = Jm«r). Let 
Gmi^,z) = /o°° r ij,ri{r, z) dr 

and HmiG-) = /o°° i^m{r,z) Jmi^r) dr 

for 0 < r < 1 

In that case, Gm-, Hm are the solutions of the partial differential equations, 

- e?Gn,{i,z) = Oand(^ - e) = 0 

In obtaining this solution, it has been assumed that 
(r"|^)r=i = 0 for n = 0, 1, 2 and 3 
and (r”‘|^)r=i = 0 for m = 0 and 1. 
where the radius of the cylinder is chosen as unity. 

As d and ip are unknown arbitrary functions, a suitability chosen function can be 
assumed to satisfy these conditions, e.g. 

(p = A{9) {[ogrY ox Ip = B{9) (log r)^ would satisfy. 

The general solutions are 

Gm{.^-,z) = {Am + BmZ)e^^ + {Cm + DmZ)e~^^ 

Hm{^,z) = + Fme~^^ 

As for a long wire 2 -> oo, the quantities must vanish, hence, the constants Am, Bm and 
Em are not admissible. 

Finally, cpm and ipm may be obtained by means of Hankel inversion theorem. However, 
it should be clear that the inversion is not automatic. The complex value of ^ with the 
substitution of C = will change the situation, in addition, because of the finite nature 
of integral 0 < r < 1 the inversion will be obtained in finite range only. The final result 
will be obtained from an integral equation, through Stieltjes moment integral theorem. 
The boundary' conditions are; at r = 1, cr^r = 0 and cj^z = 0 
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After simplifying one gets, 

Dm = 0, Frn = aiici also ^ Hiust Satisfy 

[J,nMO - ■/m-l(s")] = 0 

For example, wlien in = 1. Di = 0, Fi = , the roots of Ji ((^) = 0, ^ ^ 3 ,- 

so on are the roots Thus, Sriedden’s theorem [26] can be applied taking sum of all the 
roots in the finite range inversion formula Lastly, one would get the displacements as, 

1 r 

“ = O ^m=o,i,2, / Fm 3Jm-l{^T))d^ COS mO 

Z Uo 

1 r 

^ = O ^^=0.1.2, / Frn 3Jm-i{^r))d^ sm md 

Z L^O 

' roo 

w - -E“^oi,2,. / Fm ^^{Jm{^r) d^ COS md 

uo J 

The stresses can be summerized as 

COS m9 

’ ’ JO [ zr Z T 

roo 

= 2G E^,^o 12, . / Fm 6“^"^ Jm(^r) d^ COS md 
’ ’ JO 

■ roo 

Gsz = 2GE“=oi,2, . / Fm Jm-i{^r) + d^ sin md 

Uo J 

' rOO 

azr = 2GE~^oi2 / Fm Jm-i{^r) d^ cos 

JO [ zr It 

sinmd 

On = 0 plane, under the action of a normal pressure opening po at the origin point 
the deflection is equal to Wq. For Mode I the sress intensity factor can be obtained from 
the following mixed boundary conditions: 

COS md^ 

= Wq 1 — in E space (i.e. within the semi-elliptic region ) 

[ m a^j 

= 0 in z - E space (i.e. in the rest of the cross-section) 
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and 


<7,.,. = -2GE~,„,,,,, 


•4n-l(4^) 


7n + l r f 3m - I ^ 

fr> ^"'”‘V _/;7D^ 1 / fri X 


4(5,) 


v2H^ 


FmiO 


= — po in E space, (i.e. within the semi-elliptic region) 

= unknown m z-E space, (i.e. within the rest of the wire cross section) 


where, the co-ordinate systems in (r^d^) are explained in Fig 3(c) 

For Mode II deformation, the stress intensity factor can be calculated from the mixed 
boundary condition as on =0 plane, proceeding in the similar manner. 


u; = 


— — E°° / 


1 fF 

2 ^ 


§ r ( ) 

2 


Fm{0 


cosmd^ 


= £/o 


i 


1 _ ^ 
C2 0 ? 

0 in Z - E space. 


in E space 


and 


CTyl rl 


2G s: 


«l,3,5, I 






■m' 2(;fe) 


\/2fl 


3 m — 1 




j (IiL\ 

rl \ “'m-U fFTT^i 




Fm{F) ^ d^cosmd^ 


= — poin E space 
= unknown in Z-E space. 


where, F = \/2R — Fi, po is the shear stress acting on the crack to create Mode II 
S I.F. and Uq is the sliding displacement at the origin, of the upper face w.r.t. the lower 
face. 

Let us first solve Mode I S.I F. problem to evaluate Kj. Thus, solving the integral 
equation F’m('f) is to be determined to obtain a^izi on 9^ = 0, at the tip of the 

semi-ellipse. Then a similar procedure can be followed for Mode II loading. 
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Expanding both sides on the powers of i'l and finally introducing, 
fo°° = /u, with n = 0, 1, 2, 3, 4, . . . 

where, /(O = {0 

foi m = 1 i.e considering the first term m the senes only, with 9^ = 0. following 

Titchmarsh [19], using Stieltjes moment problem, 

FliOe = fiO = S- 0,1.2.. /o~ (-1)” dn ^ C0S(V^5) ds 


or 


m) 


= 


Lastly, expanding on both sides and equating the co-efficients of yi one can find the 
co-efficients Truncating upto five terms only, five linear algebraic equations will be 
obtained, hence, after solving, 


Wo 

Wo 

1 - 

3 

1 a2 

4i22 + 

9 


aW 

16 



■a2 

1 a‘^' 

n/2-* 

3 


aW 

J 

Wo 

1 

a" 


■ 2 ■ 


Wo 

9 a‘^ 




/o = 

h = 
h = 
h = 

The expressions for the normal stress on = 0 plane is 

(7,1,1 = -f [.4os(0,0) -b ^6’(1,0) - 5oz(0,-l) - 5oz(l,-l)] 


[^0 ^(0) 0) d" ■^0'S(1, 0) — Bq /(O, —1) — Bq d(l, —1) ■+■ Boc(0, 0) H- Bo c(l, 0)] 
where, z(0,-l) = I(0,-1) - c(0,0) and z(l,-l) = 1(1, -1) - c(l,0) 

as well as Aq = /q — ^ ^ -+• 720 x se 

and Bn = h- — dz. — Iz — — Ia — 
ailU JDo 2 12 90 56 X 144 

The symbols for the integrals are [20]: 


s(0,0) 

s(l,0) 



sin'U Jo{uTo)e du 
sinu Ji(uro)e"“''° du 
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TOO 

c(0,0) = / u~^(l — cosu) Jo{uro)e~'^’'° du 

Jo 

c(0,0) = /(0,-l) - 2(0, -1) 

/(O, — 1) = [ u~^ Jo{uro)e~'^''° du 

Jo 

2(0,— 1) = / u”' cosu Jo(uro)e““''° du 

Jq 

TOO 

c(l,0) = / u“''(l — cosu) Ji{uro)e~'^''° du 

Jo 

= - 2(1, - 1 ) 

roo 

/(O, -1) = / u-^ /o(wro)e-“’'° du 

7 0 


and 

/(I,— 1) = /o°° u”^ Ji(uro)e~“’'° du 

T'o = 

From the definition of S.I.F., the mode one Ki can be determined, at the tip of the 
ellipse, as 

Kj = limj/-,a \/2 y/y - a a^i^i = limj,_4o - v ^0 ( 1 - r^)"5 - a V2 

where, *4o = |^1 — 0.375^ + 0.282'^j 

Following Newman et al [12], after simplifying, 

Ki = po\/^v^[l-122-0.314f -0.139(f)2...]x[l - 0.354f - 0.438(1)2 ^ o.lll(f)3] 
By taking more number of terms in the series and adding higher values of m, the accuracy 
can be improved. 

Proceeding in the same manner, for m = 1 and 6^ = 0, the stress intensity factor 
for Mode II can be determined, after expanding, and then equating the coefficients of 
yi, one obtains similarly. 



/o 

ii 

h 

h 



h = 


11 



h 


0 


/e = 0 

higher powers of have been neglected 
Thus, finally, for a Mode II shear crack, 

Mo = Uo[l - 0.138^ - 0 381(|)''' 

After simplifying, and taking limits, the stress intensity factor for a Mode II shear 
crack, at the tip of the semi-ellipse becomes, 

K[[ = go \/^\/a[l.l22-0.211(f)2 - 0.094(a)^]x[l - 0.354| - 0.201(|)2 + 0.0267(^ 
Accuracy can be improved taking higher order terms, po and go are the shear stresses 
acting on the semi-elliptic part-through crack. 

An inclined crack under biaxial loading, the stresses near the vicinity of the crack 
tip are given by [15], 

^yy = ife I [l + sm f sin f ] + ^ sin f cos f cos f 

<Jxx = cos| 1 — sin I sin Y — sin | j ^2 -f- cos |cos y] + <^(I — cos 2 a 
sinfcosfcosf 4 - ^ cosf [l - sinfsinf] 
where, (r, 9) are usual polar co-ordinates with respect to the crack tip and crack plane, a 
is the crack orientation with respect to the vertical loading axis and K being the biaxial 
loading factor. Fig. 4(a) shows constant shear stress contour lines observed from these 
equations [ 11 ]. 

In our problem, radius of the wire = 1.0 mm, a = 0.11 mm, ^ = 0.11, ^ = 0.25 
from the experimental data. Thus, if one employs the angle of initial crack extension 
under biaxial load conditions then from [15] 

{^§^)e=eo = 0 OT 

l.7b{B\{a,K) [3cos^o ~ 1] ~ ^B 2 {a, K) sixig^j} -t- 4 ^ 3 ( 0 !, AT) sin^ cos^o = 0 
where, (a, AT) = —3(1 — A") sinacoso; ^ 2 ( 0 ;, AT) = AT 4- (1 — ii")sin^a 
and S 3 (a. A) = (1 — A) cos 2a 

.■\ccording to Verpoest [1], ^ = 0.5 on 45° inclined plane, then from Morh’s circle, 
biaxiality factor A = | taking 1 / = after solving the above equation, one gets 
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^0 = 54.9°, which means the crack-extension is along 9.9° inclined with respect to the 
axis of the wire, which does not agree with the experiments 

Following Irwin [7], Sih’[8] and Smith[9]. if the strain energy density criterion is 
employed, then. 


r — S — CliiKf -|- 1 i! 


II 


where, ^ is the strain energy density per unit volume, which is usually expressed by 
TT = + CTyO, 4- a-a,,) -b + r ^ J 

“11 = - 4i/ - cos^)(l 4- COS0)] 

ai2 = Y^2sin0[cos^ — (1 — 2u)] 

022 = lie ~ cos^) 4- (1 4- cos^)(3cosi9 - 1)] 

In order to have a crack extension Sih’s [16] criterion 
= 0 IS to be satisfied. 

For 1 / = 4, this equation simplifies as 

^ - 3%^ cot 2^0 - I = 0 

Solution of this equation gives, angle 9c at which the crack extension should occur. 
Taking ^ = 0.5, ^ = 0 11 and ^ = 0.25, solving the equation, 9o = 45° + 3.3° 
oriented with respect to the axis of the wire. Fig. 4(b) shows this orientation angle, 
calculated from elasticity theory, which does not agree with experiments. 

It will be demonstrated experimentally that elasticity theory alone is inadequate to 
describe stable crack extension in textured materials. Thus, there is a need to calculate 
the plastic stresses within the yielded medium in the presence of anisotropy. 
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1.4 CALCULATION OF PLASTIC STRESSES 


This part of the paper mainly describes the method of anisotropic plastic stress calcu- 
lations by slip line theory iii a strain-hardening solid in the presence of non-homogeneous 
plastic flow In Cartesian co-oidinate, the equilibrium equations are, 

^ + %" = 0and%- + |f = 0 
Hydrostatic pressure: p = —^{cXx + ay) 

Pure shear stress in yielding K = 


\{ax - ay)^ + 


The direction of maximum shear stress 6 = \ tan ^ — 

Hence, a^ = —p — K sin 2^, Oy = — p -f- K sin 20, and r^y 
The continuity condition is H- |^ = 0 
The directions of principal stress and strain must coincide, 


K cos 20 


2ti 


du I 


= Itf 


dv 

dx 


~1I 

dy 


tan 20 


Let a and (3 are the well knon slip lines for a rigid-perfectly plastic solid. For strain 
hardening plastic solid, the third characteristics refers to the space differentiation along 
the direction of particle flow; the angle 0, measured between the stream line tangent and 
the direction of the maximum shear strain rate direction (a). The equations referred 
to the characteristics for steady plane flow of a strain-hardening plastic material, the 
equilibrium equations are [21]; 


4. 2ir ^ ^ 

da da dp 

^ _ 2K^ - — 

dp dp da 


= 0 


= 0 


Geinngers’s mass continuity equations are; 


du 

da 

dv 


d(f> 

^da 


+ u 


dp 


0 

0 


and the hardening equation is; | Ag — sin20.AiT = 0 

where, q = \/u^ = the speed of flow and 0 is the direction of flow relative to the 

a slip line direction, 
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i.' = tan - 0 

The hardening modulus ^ ^ ^ ^ 

These equations constitute a partially- hyperbolic quasi-linear system satisfying the 
isotropic hardening rule, which also implies that, the stress and strain-rate tensors are 
co-axial. The incompressibility condition 
+ £;/ = 0 IS automatically satisfied 
The effective strain increment is given by 

= I + del + ^ 

^ , dq 

esc 2ip . — 

Q 

Thus. /i = I of the slope of uniaxial stress-strain curve. The effective strain may be 
obtained by integrating along the stream line; 

e = ^ J CSC 2i/i dg 
which is known as Stieltjes integral. 

Hence, increment of effective strain in terms of the directional derivative of q along the 
streamline can be calculated from the above equations. 

Now, in reference to the local axes along the streamline tangent and normal direction, 
this strain increment can be resolved into direct and shear components, 
d^t ^ ^^,des = q~^ cot 20; dq 

Finally et = f cot 2^ dq 

The directional derivatives appearing in the equilibrium equations can be computed by 
numencal differentiation with suitable approximations to the arc lengths. The equilib- 
rium equations are now integrated along the slip lines to give the mean pressure p, apart 
from various other additive constants. Effectively, there will be two matrices, viz., •pi 
and p 2 - 

Therefore, the equations of rigid work hardening plasticity from a hyperbolic quasi- 
linear system with three families of characteristics. The third characteristic originate 
from the streamline. 


2 . ,2 
—7= Cx CSC 20 dt = —7= 

Vs Vs 
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Thus, micro-gnd distortion study gives some clue how the plastic shear and di- 
latational strains are distributed. Chequerboard microgrid patterns, printed on the 
polished surface of the specimen, at the tip of an yielding crack, helps to study the 
non-homogeneous plastic flow field. 

A simple numerical scheme has been drawn calculate the stresses, (a) Smothening of 
the flow pattern, (b) Computation of strain rates and the streamlines, (c) Construction 
of slip line mesh, and calculation of velocity field, (d) Computation of nodal points, and 
finally, the calculation of stresses. 

The purpose of carrying out computation in this manner is not only to calculate the 
stresses, but also to identify the mismatch field m discontinuous deformation pattern 
[21] . The radial and tangential stresses are calculated from; 

= K cos 21 /) and Opp = p ■¥ K cos 2ip 

Suppose that, a rigid-perfectly plastic solution is known for a given configuration, 
then the magnitude of the velocity jump across certain boundary will also be known, 
and hence, the effective shear strain jump [7] is determined (where 7 = ;^ ). In the 
case of hardening plastic flow, taking h as small, the change in 7 across the strip is 
A7 = [7] approximately. Hence, one gets, 


aa: 

dK 

da 



where, b denotes the strip width and all quantities are to be interpreted as average of 
strip. 

The main objective is to obtain a sufficient degree of smoothning to enable the 
subsequent computation of the slip line directions, as the directions of maximum shear 
strain rate, which in turn, can be used to construct a slip line field, for further analysis 
of the plastic flow. Accordingly, a maximum step size of about 0.47 /r was used, because 
this was the pitch line gape of the micro-grids. 

The stress calculation involves estimation of the mean pressure p and the shear stress 
K at the slip line nodes, and their integration round a suitable path to obtain the stresses. 
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The various stages of computation are; (a) interpolation of e. and onto slip line 

nodes, (b) determination of K and the directional derivatives and ~ on the basis 

of an experimental stress strain curve, (c) determination of the velocity field, and (d) 
integration of equilibrium equations 
The strains are calculated as follows 


da 

dt 

dp 


= Cx cos(f) + Cy sm4> 

= —ex sin0 + ey cosp 


where these directional derivatives satisfy the boundary conditions, and ey are the 
normal strain components along the direction of x and y. equilibrium equations are then 
integrated along the slip lines to obtain p, apart from an additional constant. Hence, 
p = ~2!Kd4, + !%ds 

where this integration was carried out on an a -slip line. 

Thus, one would get the matrix P 2 of nodal hydrostatic pressure, similarly, a matrix 
Pi is computed by integrating along /? lines. Discontinuities divide the entire region into 
subdomains. By patching together local solutions for all subdomains, a global solution to 
a given boundary value problem can be constructed. In hardening flow, the slip bands do 
not propagate jumps in the velocity, the passage from one Hencky-prandtl net to another 
hardening slip line field must then include replacement of strong discontinuity', because, 
a real strain-hardening plastic material does not have a sharp stress discontinuity like 
a rigid ideally plastic solid. During stress calculations, the sharp stress discontinuity is 
replaced by a gradually increasing stress elevation smoothly changing over the region in 
a narrow band. The best possible way to do such opening a fold is to apply Whitney’s 
theorem. Fitting a function of -f ep, it would be possible to open the fold, provided 
a suitable value is taken for e with proper sign. Local corner solutions with associated 
slip line field has been constructed from experimental observations . 

k completely graphical procedure for constructing an approximate plane strain slip 
line field for anisotropic material is followed. The slip lines of anisotropic material under 
plane strain conditions are the trajectories of maximum shear strain-rate. Fig. 5(a) 
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shows that the a -slip line is inclined to the x-axis (in anticlockwise sense) at an angle 
cp. Then the outward normal to the yield locus is inclined at an angle 20 to the axis. 
Rice [22] has showen that, when referred to curvilinear a, 3 lines, Hencky’s equilibrium 
equations can be integrated to give 
p -f s = a constant along a -line 
p - s = another constant along /3 -line 

where, p is the mean compressive pressure and S is the arc length around the yield locus, 
measured in the anticlockwise sense, as shown in fig 5(b) 

In Figs 5(c) and 5(d) these equations have been applied around the curvilinear 
quadrilateral ABCD. After eliminating the pressure, 

Sb - Sa = Sd — Sc and Sc - Sa = Sd - Sb 
Thus the change in the arc length between two curves of one family, at their points 
of intersection with a curve of the other, is the same for each intersecting curve of the 
second family. The details of the method of graphical construction can be obtained from 
Rice [22], Chitkara and Collins [23], Johnson, Malherbe and Ventor [24]. 

The plastic zone size can be calculated from Hill’s [25] anisotropic criterion; 

In general, 

ai(cri - ay^ -b a2(o-y - a^)^ -b 03(0-^ - -f a^r^y -b + aer^^ -b 1 = 0 

For steel wire, fii cr^ + /j ,2 CTy + fiz cTxCTy + + 1 = 0 

Where the co-efficients are given by Mols [6] following Hill. 

Therefore, inside the plastic zone, the slip line field theory proposed by Oxley [21], 
Hill [25] and Rice [22] for anisotropic material has been applied. 

In the case of rigid-perfectly plastic isotropic material, the method of slip line can be 
used to determine the stresses. Henckey’s equilibrium equations and Geiringer’s velocity 
equations are used to find the direction of characteristics. In isotropic material the 
direction of maximum shear strain rate coincide with the slip lines, but for anisotropic 
material it is not so in general. The inclinations ^ of the slip lines with respect to the 
x-axis of anisotropy can be shown to be the roots of equation; 

((7x — (7y){dx^ — dy^) + 4(1 — c)rxj, dx dy = 0 
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And following the usual notation, reflecting to the slip lines as the a and J lines, the 
equilibrium equations can be formulated as 
^ — constant on a -line 

^ -1- g{(t>) = constant on /? -line 

In the above equations, p = + crj) = the mean compressive stress, and 

/i\ “ Sin 20 cos 20 , 

9{(t>) = 77^ + /o^(l -c sin^2<?!» 2 dcf) 

(1 -c sin^20)2 

when c = 0, this equation reduces to isotropic case following Henckey's relation 

The velocity relation along the slip lines can be shown to be as follows; 
du — vd4> = 0 along an a -line 
dv — ud(p = 0 along a p -line 

where, u and v are velocity components along the a and P lines. These are the same 
as Germger’s equations for an isotropic material. Fig. 5(e) describes the method of 
constructing slip lines in an anisotropic plastic solid. 
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1.5 EXPERIMENTAL 


Experimental investigation on crack propagation is carried out on 0.6 percent carbon 
steel wire which was cold drawn with a moderate reduction to get a finished diameter 
of 2 mm On the contrary, crack tip plastic flow is studied on a 4 mm thickness plate, 
made of 7075-T6 aluminium alloy, (which normally contains 5.5% Zn in dispersed particle 
form). Aluminium plates are heavily cold rolled along the direction of pin-loading, (see 
fig 2(b)), so that, the grain flow orientations are always in parallel to the direction of 
tension being applied. This steel wire has an yield strength of 750 Mpa with ultimate 
strength as 900 Mpa, whereas, the percentage elongation at fracture is in the order 
of 12% The major objective of this investigation is to understand the role of plastic 
anisotropy on the mechanism of micro-fracture initiation within the plastically yielded 
crack tip region In order to study this problem, steel wires are loaded in a combined 
mode of tension and torsion. A fine spark-erosion cut was given on the surface of the 
wire to initiate a fatigue crack. Fatigueing was carried out in the presence of tension 
applying alternating torque in a 10 Ton Instron machine at a low cycle range (0.5 to 0.8 
c/s). During fatigue cracking, A/f, as well as, axial tension, was maintained in the same 
way as suggested by Verpoest, Aernoudt et al [3 — 5]. After 20,000 to 30,000 cycles of 
fatiguing, the wires pulled in tension. 

Fig. 6 shows grain flow structure of the original steel wire before loading. Grains are 
seen to be elongated along the along the axial direction of the wire. Figs. 7 to 11 show 
scanning Electron Microscope photographs for fractured surfaces of the wires, which 
demonstrate the presence of a fatigue crack emanating from the notch. After careful 
study of several samples, it is concluded that, a semi-elliptic fatigue crack surface is 
initiating from the wire edge. It is also observed that, this crack is trying to propagate 
through the wire along 45° plane in a helical manner. Figs. 2, 3, and 11 agree with 
the experimental observations of Verpoest, Aernoudt et al [3 — 5]. Thus, it is felt that 
there is a need to study the crack propagation behaviour of anisotropic steel wires under 
a combined mode of tension and torsion loading. Figs. 7 to 11 show that the fatigue 
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crack first starts in pure mode I; then it grows along 45° inclined plane through the wire 
in a helical manner, under the action of combined mode I and II, and finally, the crack 
propagates abruptly along the axial direction of the wire causing shear fracture. 

Plastic stress calculations are based on the experimental study on flow pattern in 
7075-T6 aluminium alloy. Fig. 12 explains the role of slip lines ahead of a notch under a 
pure bending action. This is taken from Cottrell and Knott [26], Green and Hundy[27], 
Hill [25] and Johnson, Sowerby et al[28]. Fig. 13 explains the formation of a plastic 
zone at crack tip in 7075-T6 aluminium plate. This shape (in Fig. 13) is observed by 
a Differential Interference Contrast microscope on the surface of a pin-loaded tension 
specimen (see Fig. 2(b)) at 8.34 KN static load. To study slip line formation within the 
plastic zone SEM observation was found to be very suitable. Figs. 14 to 16 confirms 
the presence of shear lines within yielded zone, hence- forth, will be called as slip lines. 
These shear lines are carefully drawn to understand the slip systems within the plastic 
enclave ahead of an yielding crack tip. Fig. 17 shows the asymmetry in the plastic zone 
shape due to anisotropy, and also describes the nature of slip line formation regions. 
This kind of slip line plot is entirely based on the experimental observations on 7075- 
T6 aluminium alloy The next step is to see how Sih’s model [29] operates within the 
plastically yielded enclave. If one assumes that a crack tip model can be simplified as 
an incompressible plastic medium, fitted into a cut-out empty space in an elastic plate, 
then shrink-fit type on-load plastic stress within the enclave can be studied by a simple 
photoelastic model. Fig. 18 shows an incompressible polyurethene rubber (i/ = 0.5) 
fitted into a cut-out space inside a CR-39 photoelastic plate (u = 0.34). Not only both 
the materials were made stress free initially, but also, after joining, they are again made 
stress free. The bonding was perfect, having no residual or shrink-fit stresses present 
due to araldite adhesive joining. In Fig. 18, approximate dimensions are: maximum 
horizontal width of the rubber insert= 25.4 mm, minimum depth of the rubber insert 
under the vertical slit = 4.0 mm. Depth of the CR-39 plate = W = 55 mm, span = 4 W, 
thickness of CR-39 plate = 6 mm, vertical slit length of the cut = 10 mm, load applied 
= 99 N under three point bending. Fig. 19 explains the contour lines for constant 
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shear stress obtained from isochromatic fringes. The geometry of the poly- eurethene 
rubber insert has been chosen close to the plastic zone shape normally seen in a metal. 
Photoelastic experiment would certainly help to understand the differential mismatch 
stress (e g shrink fit stress type), created within the enclave Fig. 19 clarifies that there 
are two distinct regions within the enclave. ” Region A” is mainly due to the differential 
compression on the crack extension plane arising out from the boundary contraction 
below the slit. Point B m Fig. 19 gives a pure shear region, which is mainly due to 
the combined action of horizontal sidewise and vertical compression. This type of strain 
field is due to tension-compression changing region, commonly seen at the center of a 
bending bar. Thus. "Region B” of Fig. 19 is due to a pure distortion, which agree with 
Sih’s model [29j. 

Figs. 20 to 23 show printed (or marked) grid patterns on crack tip surface in 7075-T6 
aluminium alloy specimens. Firstly, fatigue crack was grown at Ki = 9Mpa-\/m in a 
pm- loaded tension specimen of 7075-T6 aluminium alloy (see Fig. 2(b)). Secondly, the 
specimen was pulled in tension to see the crack opening as well as plastic flow. Finally, 
grid distortion was observed under a differential interference contrast microscope as well 
as on scanning electron microscope. Photo offset grid patterns are made at Agfa-Ghobert 
factory at Mortsell, Belgium with grid size as 24/i,which was found to be inadequate to 
understand the flow field within the yielded medium. Micro- hardness indentation marks 
were made on the surface of aluminium alloy plate with a size of 6.5/i as diagonal length 
using 136° diamond p\Tamid indentor in Leitz microscope. Figs. 21 to 23 show that 
6 op grid was also inadequate. Although, in the previous work, laser beam marking 
'was employed with 0.6/i pitch distance [30], this time transmission electron microscope 
calibration piece supplied by polaron Co., England was directly put on the crack tip 
to measure plastic strain. These micro-grids are on a thin film carbon replica which is 
mounted on a copper foil mesh as a support, such calibration grids are fixed onto the 
crack tip by Silastic rubber (silicon type) bonding glue. A very thin layer of glue was 
put on the surface of the specimen, so that, a large plastic strain can be accommodated 
on the replica. Fig. 24 shows SEM photograph for such polaron micro-grids, having a 
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dimension of 0 47/,i as pitch distance between the two consecutive lines, and 0.06^ as 
line thickness. Fig 24 is taken when no load was on the specimen. Fig. 25 shows 
micro-grid distortion at crack tip due to 8 34 KN tension load on the sample. Radial 
flow pattern below the tip is confirmed from Fig 25, which means, metal extrudes out 
from the root of the crack. Fig 26 shows a micro-grid distortion pattern on 45° plane, 
giving a shear deformation Shear strain distribution is calculated from these figures. 
The crack was further loaded to 11 KN and again micro-grid distortion was recorded 
at different regions. Figs. 27 and 28 show such grid distortions, which would certainly 
help to calculate strains on the plastically yielded zone. Fig. 29 shows more radial 
plastic flow at higher load, before the onset of fracture. Thus, from these photographs, 
it is clear that, extrusion and intrusion of metal play the most important role in micro- 
crack nucleation. Fig. 30 gives strain localisation within the plastically yielded medium 
ahead of a crack tip. Such non-homogeneous plas tic flow is seen to occur at a higher 
load, before the onset of crack extension. Shear bands are seen to form within the plastic 
zone. Local strain concentration is also visible from Fig.30. Intrusion-Extrusion of metal 
must be responsible for such type complex folding deformation. Strain distribution has 
been measured within the shear band. Such flow-localization creates stress concentration 
ahead of an yielding crack tip [30 — 33]. Fig. 31 has been taken at 18000X magnification 
at crack tip. Hence, it is possible to measure all the three components of plastic strain 
at a point on the surface of a crack tip. This study has certainly helped to understand 
non-homogeneous plastic flow in an anisotropic material. 

1.6 DISCUSSION 

Our experimental investigation, described in Figs. 7 to 31, has one thing in common 
that always there will be a shear strain localisation within the plastic enclave, before 
the onset of abrupt crack extension. 

Photoelasticity experiments (mentioned in Figs. 18 and 19) as well as micro-grid 
distortion recorded on aluminium alloys (see Figs. 25 to 31) have certainly helped to 
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understand the complex strain field inside a plastic zone, particularly when, there is 
d mode transition in crack propagation for a textured material Experimental investi- 
gations on crack propa gation on 0 6/4 carbon steel wire, and crack tip strain measure 
meats m aluminium alloy also explain the role of grain flow direction on micro-fracture 
initiation in steel Figs. 7 to 11 further elaborate, the transition in crack propagation 
modes, for drawn steel wires, which were fractured after giving a number of cyclic loading 
in a combined mode of tension-torsion. 

Recently, attention has been focussed on the propagation of an inclined crack un- 
der combined mode of loading [10, 11, 13, 15andl6], but plasticity aspect has not been 
considered properly, especially, in the background of new development on stable crack 
extension theory. This work may help to under stand sub-critical crack growth for an in- 
clined crack in a textured metal. Although, Irwin-Sih-Smith approach [7 — 9] has helped 
to find conditions favourable for crack-extension in elastic materials (like polymer, rub- 
ber or glass), unfortunately, such calculations may not be aplicable in metals, especially 
without any previous knowledge on plastic stresses or anisotropic grain flow conditions. 
Stress concentration due to strain- localisation makes the situation more intresting, in 
the presence of anisotropy. Non-homogeneous deformation and reverse plastic flow with 
shear localisation within the plastic zone are seen to occur. Figs. 7 to 11 demonstrate 
abrupt crack-path instability for one inclined crack in a steel wire, because of anisotropy. 

For a cold drawn steel wire, a crack initiates purely in mode I (see Figs. 10 and 11) 
as expected, and then, propagates in a combined mode with Mode I and mode II. Fig. 
11 shows that the crack- growth is seen to occur along 45° plane in a helical manner 
under combined mode [1]. For ths 45° inclined semi-elliptic crack, it is emphasized that 
the shape of the plastic zone will always be asymmetric. This is mainly due to the 
mixed mode loading on the crack and plastic anisotropy [6, 11, 15] (see Figs. 4(b) and 
11). Next, two distinct conditions are to be satisfied: firstly, one of the rabbit shear- 
edr’ of plastic enclave should be parallel to the axis of the wire, and secondly, a critical 
shear stress concentration would be needed within the plastic zone [30 — 34]. Thus, at 
some stage of the growth, this inclined semi-elliptic crack suddenly bifurcates along the 
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axial direction of the wire This is mainly due to the sudden loss of resistance to shear 
deformation parallel to the axis of the wire, and happens to be the anisotropic gram 
orientation direction of the texture Such unstable crack growth can easily be observed, 
if the lesolved components of stresses are chosen in such a manner, so that, on 45° 
inclined plane, the applied shear stress is half of the normal stress being applied [1], i.e 
^0 = 0 5po- 

Our calculation based on elasticity theory suggests that a 45° inclined semi-elliptic 
crack should extend along a direction of 8.3° with respect to the aocis of the wire. How- 
ever, experimental observations (see Figs. 7 to 11) on steel wires do not agree with this 
The discrepancy is due to the fact that , neither {Tmax)0=:8o) is a critical shear stress con- 
centration at crack tip for a metal, nor crack would ever nucleate at the elastic-plastic 
boundary. The objective of this work is to prove that , in order to have a crack extension, 
the local maximum shear stress vector within the yielded zone must coincide with the 
direction of anisotropic texture 

In a drawn steel wire, not only one should know, the asymmetric shape, size and 
orientation of the plastic zone for an inclined crack, but also, it is equally important 
to determine the maximum shear stress concentration and its direction with respect to 
anisotropy. A simplified approach based on the findings on aluminium alloy has been 
presented here, because, it was not possible to put micro-grid on a steel wire. Figs. 27 to 
31 give more non-homogeneous plastic flow due to anisotropy, compared to our previous 
investigation on isotropic 7075-T651 aluminium alloy [30]. Figs. 32 presents the final 
result for for the shear stress concentration on 45° plane, ahead of an yielding crack tip 
in 7075-T6 aluminium. A possible slip line field is conceived in Fig. 33 [30]. 

Thus, for an inclined crack in a wire under combined mode of tension torsion abrupt 
change in the direction of crack propagation path is due to the followdng reasons: firstly, 
one of the “rabbit ears” of the plastic zone must orient itself along the axial direction 
of the wire, secondly, the resolved maximum shear stress concentration on 45° plane, 
within the plastic zone, must reach at the critical value, and lastly, at some suitable 
point inside the material, during a favourable stage of propagation, the magnitude as 
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well as the direction of this maximum shear stress must coincide with the minimum 
resistance to shear deformation offered by the anisotropic wire 
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Fig. 1(a) Cold-drawn cylindrical steel wire of 2 mm dia under combined mode of loading in 
tension and torsion. Grain flow direction is parallel to the wire axis A spark erosion 0 1 mm 
notch cut is given at the middle length of the wire, as fatigue crack staiter 
Fig. 1(1)) Olihqiie necking in a strip under tension [25], Direction of gram flow is not parallel 
to the diiectiou of tension loading Mateiial, 7()75-T6 aluminium alloy plate 




Fi{>, 2(1)) Diniciisions for a typical pin-loading type tension specimen as per ASTM Maleiial, 
707o- 1 6 alnininiiiin alloy. Giain How dnection is jiarallel to the tension loading axis 
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Fig, 4(a) A 60 inclined crack is in a tension plate. Constant shear stress contour lines <ue 
calculated from References [11, 15, 16] to obtain photo elastic data for different fringe values 
Fig 4(b) Plastic zone formation is explained at the tip of an inclined crack in a circular wiie 
T,na,. at 0 = 0i) is also shown Irwin-Sih elastic theory suggests that the angle of crack extension 




Fig 5(a) Relation between the slope of slip-lines and normal on yield locus [231. 
Fig 5(bj Arc length measured on the yield locus. 

Fig 5(c) Slip line construction by graphical method. 

Fig. 5(d) \ariation of angle <p on the yield locus at different net points. 

Fig 5(e) Final modified shape of the slip net along a and 0 directions. 
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Fig. 6 Gram flowstructure for 0.6 C steel wire, obtained after cold drawn through the 
HSS die. Grams are seen to be elongated along the axis of the wire. 
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Fig 7 A typical helical fractured surface of a wire after fatigue cracking. Scanning 
Electron Microscopic observation at low magnification. 
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Fig 9 SEM photograph taken at an angle showing different fracture surfaces for a 2 mm 
steel wire. 
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Fig 10 SEM photograph showing the initiation of a surface fatigue crack from spark 
notch cut. 
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Fig 11 Same SEM photo of Fig 10 on the right hand side, but a small region is 
magnified XlO times more than the right hand side, showing fatigue crack on the L.H.S. 
figure. A semi- elliptic crack is foiming on elliptic plane inside the wire [5] . conditions 
tbi the crack growth is taken from Ref [5], but grown in a low cycle fatigue machine. 
^Paxiai = 0 IFfracturc , AT = 0.85 Tfracture, ^ = 28,000 cyclcs, where T is torque 


and P is the axial load on the wire 


Plastic zone on 0 plane 

Negative velocity Positive velocity 

vector L vector 



Fig 12 Slip line field from a deep groove V-notch under pure bending action on the 
beam. Experimental observation by Fry’s reagent on Nitrogenated steel is also shown 


on a photograph [25 - 28]. 
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Fig 13 Plastic zone measured on the surface of 7075-T6 aluminium fatigue cracked beam. 
Photograph was taken in a Differential Interference Contrast microscope. Crack size = 
3 41 mm. Pin loaded tension specimen. Load applied = 8.34 KN. Pre-cracking for fa- 
tiguing A AC = 9Mpa with R = 1 AACc = 27 M pa -^/rn 
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Fig. 14 Slip systems observed at crack tip by SEM study, conditions are same as Fig. 
13 Material; 7075-T6 aluminium. 
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Fig. 15 Same slip systems observed at crack tip of Fig. 14, but a higher magnification 
(7075-T6) 
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Fig. 18 Photo-elastic isochromatic fringe pattern is shown inside polyurethene rubber 
insert, fitted inside a CR-39 photoeiastic matrix. Differential mismatch loading is clear 
inside the rubber insert. Load = 99 N, 3 point bending loading, CR-39 specimen sizes 
are W = 55 mm, S = 4W, B = 6 mm, L = 250 mm, slit size = 10 mm minimum insert 
depth under the slit = 4 mm and the maximum width of the rubber insert = 25.4 mm. 
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Fig. 19 Constant shear stress lines inside the rubber insert due to differential mismatch bound- 
ary movement. ’Region A’ demonstrate the presence of a differential compression, whereas. 
■Region B’ confirms a pure shear field due to bending and side- wise compression. This helps to 
understand the slip lines inside plastic enclave. This agree with Sih's model[29] in metals. 
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Fig. 20 grid distortion study at crack tip, before the onset of fracture. Grids were made 
by photo-offset technique at Agfa-Gobert, Mortsell, Belgium. Size of each square grid = 
24 n mesh. Material; 7075-T6 aluminium alloy. Crack size = 3.69 mm. Pin-loaded ten- 
sion specimen. Load applied = 1250 KN Fatigue pre-cracking data A AT/ = 9Mpa ^ 
with R = 0.1. 
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Fig 21 Grid study on the ‘surface of a 7075-T6 aluminium alloy specimen (by Leitz 
micro-hardness tester) indentations. Grid sizes = 6 5 /j,, pitch between the consecutive 
grid centers = 10 fi. Crack size = 3.69 mm. Pin-loaded tension specimen. Load applied 
= 11 KN. Fatigue cracking at AKj = 9Mpa ^/m with R = 0.1. 
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Fig 22 SEM study of fig. 21, but taken at a higher magnification. Conditions are same 
as in Fig 21. 
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Fig 23 Same conditions as Figs. 21 and 22, but taken at a lower magnification, on the 
R.H S., and a small portion is also shown at a 2000X times higher magnification on the 
L H S. through SEM photograph. Micro-grid elongation is visible below the root of an 
yielding crack tip. (Take same data as Fig. 21). 
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Fig 24 Micro-grid photographs by SEM. Transmission Electron Microscope cross-grid 
calibration piece from polaron Co, England, is mounted on copper foils. Carbon replica 
thickness is in the order of 400 Ao- Replica square grid lines shown here, has pitch = 
0 47 fj,. Line thickness = 0.06 /j,, such T.E.M. calculation piece is glued on the surface 
of a 7075-T6 aluminium specimen near the crack tip region, by silastic (silicon) rubber 
adhesive, which can take large strain deformation. 
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Fig.25 SEM photo for a crack tip micro-grid distortion study. TEM square grid pitch 
distance = 0.47 /i Radial flow pattern is visible below the crack tip due to metal extru- 
sion Material 7075-T6 aluminium. Crack size = 3.41 mm.Pin-loaded tension specimen. 

Load applied = 8 34 KN. Fatigue pre-cracking data AKj = 9Mpa -s/m with R = 0 1 
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Fig. 2C SEM photograph of a carbon replica micro-grid distortion at a crack tip taken on 
45'^ plane with respect to the extension direction Shear strain distribution with in the 
Yielded zone has been measured on the surface of a cracked specimen. Material; 7075-T6 
aluminium alloy Crack size = 3 69 mm. Pin-loaded tension specimen. Load applied 
= 8 34 KX Fatigue pre-cracking data AKj = 9Mpa -Jm. Original square grid pitch 
distance = 0.47 p. 
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Fig 27 SEM photograph on a crack-tip with TEM calibration carbon replica square 
grids mounted. Plastic strain distribution is visible ahead of crack tip. Strains have 
been measured by grid distortion study. Mixed flow field is visible Material; 7075-T6 
aluminium alloy. Crack size = 3.41 mm. Pin-loaded tension specimen. Load applied = 
11 KX. Fatigue pre-cracking data A/\/ = 9ikfpa ^/m Measured original square grid 
pitch line distance = 0.47 p , Readings were taken near 45^ line w.r.t. crack extension 
plane. 
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Fig 28 SEM photograph taken at a higher magnification on crack tip plastic flow in 
7075-T6 aluminium alloy Strain fields are measured from such photographs taken at 
various points ahead of a tip. Crack size = 3.41 mm. Pin-loaded tension specimen. Load 
applied = 11 KN. Fatigue pre-cracking data AKi = 9Mpa y/m. Original square grid 
pitch distance = 0.47 /i. Similar observations were made at various angles at different 
points with the yielded zone. 
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Tiacing of Fig 28 


29 SI'3M pliotugiaph of caiboii replica giuls oii crack tip Radial [)la.stic How is seen 
again below the crack tip at a higher load Crack size = 3 41 mm. Pin-loaded tension 
specimen Load applied = 11 KN Fatigue pre-crackmg data AR/ = 9Mpci \/m. Orig- 
inal scpiaie grid pitch distance = 0 47 


50 


J— — — — —Mllll 

■wwmwm —mi — MI 

MBBSBMWaaai 


Fig. 30 Strain-localization with shear bands, inside a plastically yielded crack tip. SEM 
photographs taken on a surface a fatigue cracked 7075-T6 aluminium alloy specimen. 
Crack size = 3.41 mm. Pin-loaded tension specimen. Load applied = 11 KN. Fatigue 
pro-cracking data AAC = 9Mpa Original square grid pitch distance = 0.47 p. 
.\ou hoinugcneous plastic strain accumulates ahead of the crack tip.Readings were taken 
on 45° plane w.r.t. the crack. 
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Fig 31 Plastic concentration inside the yielded zone at crack tip in 7075-T6 aluminium 
alloy. Plastic strain distribution is measured by micro- grid distortion method. SEM 
photograph showing carbon replica grid distortion at crack tip. Readings were taken 
at vaiious points in different regions on the surface of the specimen. Crack size = 3.41 
mm Pin-loaded tension specimen. Load applied = 11 KN. Fatigue pre-cracking data 
AAV = 9Mpa s/m Original square grid pitch distance = 0.47 p. Non-homogeneous 
plastic flow is confirmed from the grid distortion. 
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Fi”. 32 Shf^ar stress concentration is plotted on 45^ plane with respect to the crack. Stress 
calculanon is btised on slip line field theory described in the text. Graphical construction was 
(‘lapluved. Slip lines are drawn from expermental observations. Mixed stress field is noticed 
Shear localization is due to the intrusion-extrusion of metal. Material 7075-T6 aluminium alloy. 
Te.-^t conditions are exactly the same as described in figs. 27 to 31. 
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Fig ;53 Slip lines proposed from the expermental observations on crack tip plastic flow alpha is 
fh(' wf'll known slip line due pure bending action on the notch about plane OA (i.e. 0° plane) 
0 IS another slip line due to the bending action about OB plane, beta slip line has also been 
ob.sf'ived expernientally while studying the shear deformation ahead of the crack tip. This is 
ba.sfd on the study of plastic flow at crack tip m 7075-T6 anisotropic aluminium alloy. Stress 
ccmccntrations are seen to occur at P and Q. Micro-fracture can take place in these regions. 
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Chapter 2 


WARM PRE-STRESSING 
EFFECT ON FRACTURE 
TOUGHNESS OF MILD STEEL 

2.1 INTRODUCTION 

This part of the thesis deals with the warm pre-stressing effect on crack toughness of 
mild steel. Firstly, a cracked beam is loaded at room temperature, and then, the beam 
is fractured at sub-zero temperature. Toughness is seen to increase at this temperature 
due to warm pre-stressing. The prominent theme of the thesis has been to find the 
relationship between the brittle fracture and plastic deformation. 

If one plots the uniaxial tension as a function of temperature, in ferritic steel, it rises 
rapidly with failing temperature. The classical theory of fracture of a potentially ductile 
material fracture in brittle is investigated in this chapter. The entire ductile brittle tran- 
sition behavior of mild steel is studied from looking into the effect of warm pre-stressing 
Above the embrittlement temperature cleavage fracture can still occur, if the specimen 
contains a notch or a crack, sharp enough to give the plastic constraint, superposing a 
tria-xial tension and raising the general tensile stress during plastic yielding. Thus, the 
u-arm pre-stressmg on notch embrittlement is studied at sub-serc temperature. There- 
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fore, phistic deformation is an essential part of the brittle fracture of metals. It may 
nucleate cracks of the critical size required by the Griffith fracture condition, or it may 
h(' able to propagate by a mechanism different from that of the Griffith theory. Thus, 
the strain strengthening is playing on important role in fracture. 

Once a crack is found on a cleavage plane, the only work required for brittle fracture 
IS that needed to break the interatomic bonds at the crack ends and this can usually be 
provided by the elastic energy at the specimen As the crack does not require thermal 
activation the velocity of crack propagation can be very high, about half the speed of the 
•shear waves in the material. Because of the large amount of work expended in plastic 
deformation, ductile fracture propagates at much lower velocities, and in certain cases 
fracture can be stopped. Thus, brittle and ductile fracture can be distinguished by their 
mode of propagation, even though in both cases the deformation may be of a ductile 
naturi' prior to the start of crack propagation. Many solids like mild steel, can exhibit 
both brittle and ductile fracture characteristics, depending on the mode of stressing, the 
rate of application of the load, and the temperature. 

2.2 RESULTS 

Fig. 2.3 shows the effect of fracture toughness on warm pre-stressing. Firstly, a 
cracked beam in loaded upto 70% of the fracture load at room temperature.Then the 
sample is broken at sub-zero temperature. Temperature range was from -80° C to 
-r50° C. Ductile-brittle transition behavior is seen to occur, which is plotted in the 
figure. Next, same thing was repeated at 80% of the fracture load at room temperature. 
Thus, the sample is loaded in Fig. 2.4 upto Q.SPc, when P, is the fracture load at room 
temperature(Fig. 2.5 and 2.6). Some samples were tested without warm pre-stressmg 
at room temperature, which is shown in Fig. 2.7. Fig. 2.8 shows the Kjc plot against 
warm pre-stressing effect. From 50% to 60%. there is no change in toughness. However, 
above 60% warm pre-stressmg, there is an increase in toughness.Therefore, due to warm 
pre-stressing fracture toughness of the metal is increased 
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2.3 DISCUSSION 


At room temperature, when a notched beam is loaded, the plastic deformation spreads 
at t h(' tip of the ciack The material is permanantly damaged m this region. More volume 
of material is deformed compared to the same at sub-zero temperature. Bandyopadhyay 
[.‘10 — 33] model shows the interaction mechanism, between the free notch surface and 
till' radial stress discontinuity. This means, at a larger distance from the crack tip, the 
material is damaged, and cracks, voids etc appear locally in this region. When the 
specimen is loaded at room temperature, such dajnaged material remains there, even 
if the sample is dipped into a lower temperature. The surface insttability load will be 
moK' at lower temperature due to increased interaction stress. Thus, at low temperature, 
higher crack extension force will be needed to extend the crack. Fig. 2.3 to Fig. 2.6 show 
continuous increa,se of toughness with temperature Ductile-brittle transition behaviour 
is siH'ii to occur in all the cases. Fig 2.8 shows the increase m crack-toughness due to 
warm pre-stressing effect. 





I'iu 2 2 is the uniaxial and C is the notch-constraint yield stress, S is the fiactuie 
.tu-ss (or ( hxnaRC fiactnrc The dashed cnr%e parallel to, and just al)o^e, V is the cm ve 
nf cleinaRe-liactine stiess in uniaxial tension below 11 The dotted cune shows tlu' 
lotus of ductile fiacime of notch specimens above H 
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Fig 2 4 Fracture toughness tests at different temperatures when warm pre- 
stressed at 80% of the Crack extension load” 


74 












CONCLUSIONS 


• Elasticity theory alone is inadequate to describe stable crack extension in elastic- 
plastic metals. 

• Shear strain localization within the plastic enclave is seen to occur before the onset 
of abrupt crack extension. 

• A critical shear stress concentration is needed to extend the crack in a metal 

• Part-through semi elliptic crack is seen to form in a circular wire while the crack 
propagation in helical path. 

• Calculations based on Irwin-Sih-Smith overestimates the direction of maximum 
shear stress vector by 3.3® with respect to the axis of the wire. 

• Crack toughness increases at sub-zero temperature due to warm pre-stressing. 
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